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We calculate exactly the Josephson current for c-axis coherent tunneling between two layered super- 
conductors, each with internal coherent tight-binding intra- and interlayer quasiparticle dispersions. 
Our results also apply when one or both of the superconductors is a bulk material, and include the 
usually neglected effects of surface states. For weak tunneling, our results reduce to our previous 
results derived using the tunneling Hamiltonian. Our results are also correct for strong tunneling. 
Ifowever, the c-axis tunneling expressions of Tanaka and Kashiwaya are shown to be incorrect in 
any limit. In addition, we consider the c-axis coherent critical current between two identical layered 
superconductors twisted an angle (po about the c-axis with respect to each other. Regardless of the 
order parameter symmetry, our coherent tunneling results using a tight-binding intralayer quasi- 
particle dispersion are inconsistent with the recent c-axis twist bicrystal Bi2Sr2CaCu208-i-i twist 
junction experiments of Li et al. 



I. INTRODUCTION 

One of the most interesting features of superconductiv- 
ity is the Josephson effect. This occurs when flat surfaces 
of two superconductors are brought together, forming a 
uniform junction. A supercurrent flows without a volt- 
age drop across the junction, provided that the properties 
of the superconductivity in the two superconductors are 
compatible. For conventional superconductors in which 
both the normal state and superconducting properties 
are isotropic, it doesn't matter which particular crystal 
surfaces are employed to form the junction. 

However, for anisotropic superconductors, the junction 
orientation can be very important. Even if the supercon- 
ducting order parameter (OP) is isotropic, the intrinsic 
anisotropic normal state properties of a layered supercon- 
ductor make the properties of Josephson junctions involv- 
ing one or more layered superconductors different from 
those formed from two isotropic materials. For example, 
Josephson junctions between an isotropic, conventional 
superconductor and a layered superconductor can dif- 
fer greatly, depending upon whether the junction is on 
the top or an edge of the layered superconductor. This 
is especially true if the layered superconductor has an 
anisotropic OP. In addition, the properties of Josephson 
junctions formed between two layered superconductors 
depend strongly upon the junction orientation, especially 
if the OPs are anisotropic. 

In this paper, we consider the case of a Josephson junc- 
tion formed between two layered superconductors stacked 
on top of each other along the c-axis. Our results can be 
easily modified to include the related problems in which 
one or both of the superconductors is bulk, rather than 
layered. By treating the two superconductors as layered, 
the surface states that form near the junction appear nat- 



urally in the calculation. These surface states affect the 
Josephson current results, even when the limit of tun- 
neling between two isotropic superconductors is taken. 
With one exception, these surface states have previously 
been neglected. In addition, when the quasiparticle dis- 
persions of the two superconductors are different, such 
as for different materials or identical tight-binding ma- 
terials that are rotated with respect to one another, an 
impedance mismatch occurs, reducing the coherent crit- 
ical current. 





FIG. 1. Sketch of the junction under study. Layer numbers 
and interlayer tunneling strengths are indicated, along with 
the general forms for the bare quasiparticle dispersions. 

We consider the Josephson tunneling along the c-axis 
between two distinct layered superconductors, assuming 
all of the tunneling processes can be taken to be purely 
coherent, as pictured in Fig. 1. In Section II, we solve for 
the Green's function in each layer, keeping the quasiparti- 
cle dispersions and OP symmetries fully general. The ef- 
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fects of surface states are explicitly included in the calcu- 
lations, and the various tunneling strengths are included 
exactly. In Section III, we derive the temperature (T) 
dependence of the superconducting gaps under three spe- 
cific assumptions about the quasiparticle band structures 
and three different OP symmetries in each layered super- 
conductor. These quasiparticle dispersions are chosen to 
model the experimental cases of underdoped, optimally 
doped, and overdoped Bi2Sr2CaCu208+5, (BSCCO), re- 
spectively, with respective Fermi surfaces that we de- 
note as FSl, FS2, and FS3. However, the dispersions 
for overdoped samples can also be made to fit the three- 
dimensional free-quasiparticle dispersions of conventional 
bulk superconductors. The OP symmetries chosen are 
the isotropic s-wave and the tight-binding (i^2_j^2-wave 
and "extended-s"-wave OP functions, respectively. In 
Section IV, we solve for the gaps for the symmetric case 
of two identical superconductors. Then, in Section V, 
we calculate the Josephson tunneling current across the 
two layered superconductors, assuming that the coupling 
across the junction is not stronger than the intrinsic cou- 
pling within each layered superconductor. Our final re- 
sults for the tunneling current can be used to calculate 
the c-axis Josephson critical current /c, and can be gen- 
eralized to the standard case of two bulk, conventional 
superconductors, for which the usually neglected surface 
states are explicitly included. Then, in Section VI, we 
present detailed /c('/'o) results for coherent c-axis twist 
junctions between identical coherent layered supercon- 
ductors. Finally, we discuss our results in Section VII. 

II. MODEL AND PROCEDURE FOR SOLVING IT 

We assume that two layered superconductors are 
placed with the upper (U) one on top of the lower (i), 
and that the contact between them is sufficiently strong 
that quasiparticle tunneling between them occurs. For 
simplicity, we assume each superconductor consists of 
^ 1 layers separated an equal distance s apart. We 
index the layers with n, m, where —N + 1 < n,m < N. 
In the U half space, 1 < n, m < A^, and in the L half 
space, —N -|- 1 < n, m < 0, as pictured in Fig. 1. Within 
each layer in the rj = L,U half space, the quasiparti- 
cles propagate with dispersion ^o?)(k) and have the gap 
function A^(k), where k = [kx, ky) is a two-dimensional 
wavevector. We assume is independent of layer in- 
dex. This assumption will be checked in Section III, and 
is found to be usually valid. Between adjacent layers in 
each half space, the quasiparticles tunnel with matrix ele- 
ment J,,/2. At the junction between layers and 1 where 
the L and U half spaces meet, the quasiparticles tunnel 
with matrix element J/2. Since we are not interested in 
spin-dependent effects, we assume the quasiparticles are 
spinless fermions, only taking account of the spin val- 
ues in counting the number of quasiparticles. We set 



kB = c = Ti = \. 

We shall focus on the general procedure for evaluating 
G. The details are given in the Appendix. In order to 
calculate the Josephson current across the junction be- 
tween the two half spaces, we first find the form of the 
finite temperature Greens' functions matrix G. This ma- 
trix is the product of two matrices, one of rank 2N , with 
elements indexed by the layers n, m, and the other the 
Nambu matrix of rank 2, with elements G, —G\ and 
F'^ in the usual cyclic order beginning with the upper 
left hand position. This Nambu matrix is represented by 
the Pauli matrices for i — 1,2,3, plus the rank two 
identity matrix tq. We let uj represent the Matsubara 
frequencies. 

We first begin by constructing the Green's function 
matrix Q for a bulk layered superconductor. We then 
add a perturbation with the particular form that decou- 
ples the Green's functions in each half space from each 
other. The resulting Green's functions of two sin- 
gle half spaces have parameters appropriate for each half 
space. We then couple these two half-space Green's func- 
tions together. 

The Green's function matrix Q'^ for a bulk supercon- 
ductor of type rj satisfies 

N - Co^^3 - A,,ri + JVa]^" = i, (1) 

where the layer space matrix J"^ has the elements 

■^mn = + S,n.n-1^ ■ (2) 

Terms not containing a Pauli matrix are implicitly pro- 
portional to Tq. In a bulk layered superconductor, it is 
then possible to Fourier transform this expression, as we 
did in the Appendix, with the result differing only from 
that of a bulk, three dimensional superconductor by the 
corrugated cylinder form of the quasiparticle dispersion, 
^^(k, fcz) = ^or;(k) — JjjCOs{kzs). Howcver, for a half 
space, such Fourier transformation is not permissible, 
and so we must keep the layer indices n, m explicitly. 

To construct , we first remove the tunneling matrix 
elements across the junction between the two supercon- 
ductors by adding the perturbation V^t^, with elements 

Km = - '^{^moSnl + S,nlSno), (3) 

taking account of the restrictions n, m < for rj = L 
and n, m > 1 for rj = U. The perturbation V'^ "cuts" 
the bonds connecting the two identical half spaces which 
terminate at layers and 1. We then solve for the 
elements gJJ^^ of g^ satisfying 

gv ^ gr, + g-JyVag". (4) 
Solutions for each g^ are given in the Appendix. 
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So far, we have found the expressions for two distinct 
half-space Green's functions, which are electronically un- 
coupled from each other, since no quasiparticle propaga- 
tion from one half space to the other has yet been intro- 
duced. We thus couple and together via the local 
perturbation J7 with matrix elements 

J 



Jmn — ■^('^mO'^nl + 5ml5no)- (5) 

For tunneling strength comparisons, we then let 

J2 



7 = 



JlJu 



(6) 



The exact solution to this problem of coupled half spaces 
then yields the full Green's function matrix G, with ma- 
trix elements G„m , which is given in the Appendix. 



III. GAP EQUATION 

In order to obtain the temperature dependence of the 
Josephson critical current, we need to solve a gap equa- 
tion for the temperature dependence of the gap. For this 
we make a simple assumption of a BCS-likc equation 



A,„(k,r) = |^^A,(k,k')Tr[(ri +»T2)G„„(k',c.)]. 

(7) 



w k' 



This layer-dependent gap function is in contradiction to 
the layer-independent gap function used to calculate the 
Green's functions. This must therefore be regarded as 
the first correction to the assumed constant gap, which 
determines the Green's function on the right hand side of 
this equation. Ideally, this should be a small correction, 
if our initial assumption is justifiable. 

In the Appendix, we found an expression for Goo^ts. 
Using this result, letting 



— 1 

16' 



l + -r\Eu\'\E,\^ 



(11) 



and 4)L — <t>u is the phase difference of the OPs across the 
junction. 

As a quick check of this result, we take L ^ U and 
7=1, and find that this reduces to 



Ai(k,T)=T^^A,.(k,k') 



ui k' 



„Ai(k>), 



([eoL(k')+if^L(k',C.)]2- J2 



-1/2- 



, (12) 



in agreement with the result for a homogeneous coherent 
layered superconductor. 

Note that in Eq. (|9|), the proximity effect couples the 
gap function at the interface in L to the gap function in 
U. This disappears in the limit 7 <C 1, wherein Ao(k, T) 
becomes the gap function at the surface of L: 



A(k, T)]^_surface 



rEE^^(k,k') 



Ai(k',c.)S'/(k',o.) 



(13) 



in agreement with previous results. ||^ 

The gap functions which go into the right hand side of 
are the spatially constant ("zeroth order") bulk gap 
functions, obtained from ( p^ ) for L, and an analogous 
equation for U. We use these to calculate A^(k, T) and 
A{k, T)jj,surface for T] ^ L , U , and compare the magni- 
tudes of these to the bulk results. Ideally, the difference 
in magnitudes between the bulk gaps and the interface 
gaps should be small. In the following section, we will 
investigate this for the special case of a symmetric junc- 
tion. 



eo„ + *a, + [(eo„ + ia,)2- j2]i/2 



(8) 



where ri,, = y/uj^+\A^, and defining S'^ — 5ft(S^) and 
S^' = 3(S^), the equation for the gap at the interface in 
the lower superconductor is 

ALo(k,r) 



A,(k,kM 



where 



A I" |2 t2 A 



2 



AlAu cos((/>l 



9(7) 



(9) 



(10) 



IV. GAPS FOR THE SYMMETRIC CASE 

In this section, we assume that L = U and solve for 
the bulk gap, and the gap at the interface with 7 7^ 1. 
To find the bulk gap, we fix Tc at 9 meV, and solve Eq. 
( |T^ ) for the temperature dependent gap. We assume each 
superconductor has a single OP, which for simplicity, we 
limit to s, d^2_y2, or extended- s- wave symmetry. We 
index these OPs by C = s,d,es. For each OP, we take 
the pairing interaction on the r] = L,U superconductor 
to have the form 



Ac^(k,k')-Ac,vI/^(k)*c(k') 



(14) 



where ^'s(k) = 1, ^'^(k) = cos(fcj;a) — cos(fcj^a), and 
^es{k) = {[cos{k^a) - cos{kya)]'^ + e^}^/^, where e < 1. 
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FIG. 2. Fermi surfaces studied for tunneling between two 
layered superconductors. Dashed curves: underdoped FSl. 
Solid: optimally doped FS2. Dotted: overdoped FS3. 

By fixing Tc, we determine the value of A^,;. These 
then determine the T dependence of the s-wave or d-wave 
gaps, respectively, via Eq. (|2|), and the extended-s-wave 
gap with e ^ has an identical T dependence to that of 
the d-wave gap. We take the in- plane dispersion to be 



$^o(k) ^ Jri - J,,||[cos(fc^a) + cos{kya) 
—u cos{kxa) cos{kya) — /i] , 



(15) 



where the part proportional to J^|| is chosen to approx- 
imate the in-plane dispersion relation for BSCCO. The 
value of J^ii is thus taken to be 500 meV. We choose 
three sets of parameters v and fi, which determine the 
details of the quasiparticle dispersion, and the resulting 
shape of the two-dimensional Fermi surface, pictured in 
Fig. 2. For a heavily underdoped sample, we choose 
fi = —1.3, v = 1.3, with Fermi surface denoted FSl in 
Fig. 2. For the tight-binding dispersion appropriate for 
an optimally doped sample of BSCCO, we take ly — 1.3 
and fj, = 0.6. This dispersion has the Fermi surface de- 
noted FS2 in Fig. 2. In addition, for a heavily overdoped 
sample, we choose v = and fi = 1.0, with the Fermi 
surface FS3 in Fig. 2. The only remaining free param- 
eter is the value of Jri/2, the interlayer overlap integral. 
We therefore perform our calculations for three different 
values: = 25, 50, 100 meV. In each case, = 25 meV 
gives results which closely approximate the weak hopping 
limit — > 0, as has been verified by explicitly checking 
our results with J^i — 1 meV. However, for that small an 
interlayer hopping parameter, one needs to use a much 
finer grid for the Brillouin zone integration, in order to 
obtain sufficient accuracy. 



For each OP symmetry we calculate Aj(T) from the 
symmetric gap equation. We note that Aes{T) — Ad{T) 
for Ted — Tees, which occurs for Aes = A^. In the case 
7=1, the combined junction of two layered superconduc- 
tors is just the same as a single layered superconductor, 
as long as the two halves are not twisted with respect to 
each other. However, for 7 < 1, the central junction is 
different than the intrinsic ones in each layered supercon- 
ductor. In this case, surface states can arise, and the gap 
can in principle depend upon the layer index. This 
would be particularly true in the case of two incompati- 
ble OP components, which has been considered in detail 
for a cylindrical Fermi surface previously. M 
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FIG. 3. (a) s- and d-wave superconducting gap magnitudes 
at T = as functions of 7 for FS2 with = 25, 50, 100 meV. 
(b) Plots of the T — s-wave gap magnitudes at 7 <C 1 for 
Fermi surfaces FSl and FS2, as a function of J in meV. 

For the Fermi surface FS2 and only one OP component, 
if 7 <C 1 and Jj^ = 0, then each layer would be completely 
isolated from each other, and the gap wouldn't depend 
upon the layer index, reducing to the two-dimensional 
mean-field value of independent layers. Hence, the ques- 
tion of whether the gap depends upon the layer index can 
be more easily answered by checking whether the gap de- 
pends upon the intrinsic interlayer hopping strength 
and on the impedance mismatch parameter 7. In Fig. 
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3a, we plotted As(0) and Ad(0) for = 25, 50, 100 meV, 
as functions of 7 (0 < 7 < 1), for the optimally doped 
Fermi surface FS2, evaluated with v = 1.3 and /i = 0.6, 
pictured in Fig. 2. Although there is a weak 7 depen- 
dence of Ad(0) for J,, — 100 meV, all other cases have 
essentially no significant 7 dependence, so that the gap 
values are essentially the same as for independent layers. 
The slightly greater 7 dependence of A(i(0) for — 100 
mcV apparently arises because the Fermi energy Ep is 
not much further than from the top of the quasipar- 
ticle band. In Ref. (3), it was also found that for a free 
particle dispersion within the layers that the gap at the 
surface (7 = 0) was the same as that in the bulk (7 = 1). 
Thus, our results for a finite in-plane bandwidth closely 
approximate that infinite bandwidth limit. 

We also studied As(0) as a function of at 7 ^ 1 for 
Fermi surfaces FSl and FS2, and plotted the results in 
Fig. 3b. Again, we found very little dependence to the 
s-wave gap magnitudes for 7 <C 1, so that it is generally 
safe to take the gap to be the bulk value calculated far 
from the central junction location. We also found that 
A(;(r)/A(;(0) for Q = s,d with Fermi surface FS2 were 
rather independent of J,, for = 1, 25, 50, 100 meV, each 
of the curves differing only slightly from the ordinary 
BCS curves for A(r). 



V. JOSEPHSON TUNNELING CURRENT 

The tunneling current across the junction is given by 
/ = -»eJ^(4(k)ci(k) - cl(k)co(k)) (16) 



where the Ci(k), (k) are creation and annihilation oper- 
ators for electrons with wavevector k in the i^^ layer. The 
angular brackets indicate a thermodynamic equilibrium 
average. In terms of the full space Green's functions, this 



IS 



/ = jeJTj] ^Tr[i (to + t^)(G^^ - Gqi)], (17) 



where we have suppressed the (k, uS) dependence of the 
Green's functions. 

In the Appendix, the matrices Gqi and Gio are given, 
and the trace evaluated. We find 



(18) 



where is given by Eq. (|8|) and D is given by Eq. (11). 

The equations derived above relate the supercurrent 
through the interface between two layered superconduc- 
tors to the phase difference across this interface. We can 
find the critical current of this interface by varying 



the phase difference until the maximum current is ob- 
tained. We can vary the value of 7, which we term "the 
impedance match parameter", to obtain Ic{l)- We can 
also vary T to obtain Ic{T)- Finally, we can investigate 
all of the above for s-, and d-, and extended-s-wave gaps, 
and mixtures thereof. 

To order J^, we can set D ^ 1, and // obtained from 
Eq. ( p^ ) reduces precisely to our previous result, Eq. 
(8) of Ref. (4), derived using the tunneling Hamiltonian, 
provided only that one neglects the incoherent tunnel- 
ing, and replaces J with 27o in the coherent tunneling 
part. However, this result differs greatly from that 
found by Tanaka and Kashiwaya (TK). H TK did not 
correctly derive the coherent c-axis tunneling between 
two layered superconductors, but instead attempted to 
approximate the tight-binding quasiparticle dispersion in 
the c-axis direction by treating the corrugated Fermi sur- 
face as a narrow belt around a spherical Fermi surface. 
This procedure leads to the uncontrolled approximation 
of dividing zero by zero. Thus, Eq. (100) of TK, ob- 
tained in the weak tunneling limit of their calculation, 
is not correct for coherent tunneling. Instead, it hap- 
pens (from compensating errors) to be correct for purely 
incoherent tunneling, as in the model of Ambegaokar- 
Baratoff (AB). |^ However, the subsequent Eq. (101) of 
TK is still not quite correct for any type of tunneling, and 
Eqs. (102)-(104) of TK are completely wrong for both 
incoherent and coherent tunneling. A modified version 
of Eq. (101), correct for incoherent tunneling between a 
conventional and an unconventional superconductor, was 
used appropriately in fits to c-axis tunneling between Pb 
and BSCCO. § ' 

In this notation, we can also investigate the case of 
the coherent tunneling matrix elements depending upon 
the in-plane wavevectors, J^(k^,k^) = J^(^^(k^), and 
J(kL,k!7) = J^p{\^L)v{^u)■, where (y5(k) = \ cos{kxa) - 
cos(fcj,a)|, as suggested from band structure calculations 
of Liechtenstein et al. |0] This form, along with the 
wavevector independent model, is useful for studying the 
coherent critical currents across a c-axis twist junction. 
We note that 7(ki,k[/) = 7 is independent of the k,,. 

As long as 7 < 1, the overall critical current Ic will 
be given by the above For both strong 7 = 1 

and weak 7 <i; 1 coupling with Fermi surface FS2 and 
Jn — 1,25,50,1000 meV, we calculated Ic{T), normal- 
ized to its T = value, for each of the OP symmetries 
considered. For 7^1, the s-wave Ic{T) / lc{0) curves are 
almost indistinguishable from the standard AB curve, 
as for Fig. 2 of Ref. (4). The other results are plot- 
ted in Fig. 4. Note that the d-wave and extended-s- 
wave curves are identical, but they differ substantially 
from the s-wave results. However, as shown in Fig. 4a, 
the d-wave (and extended-s-wave) Ic{T) / lc{0) curves are 
strongly dependent upon J,, , in a non-monotonic fashion, 
unlike Fig. (3) of Ref. (4). Those curves in Ref. (4) 
were evaluated using the free-particle quasiparticle dis- 
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persion within the layers. In addition, in Figs. 4b and 
4c, we plotted Ic{T)/Ic{0) for the strong coupling case 
7 = 1 for ( — s,d. respectively. In this case, the non- 
monotonicity of the d-wave (and extended-s-wave) case 
is similar to that in the 7 1 limit, but the J,, variation 
of Ic{T) / lc{0) is stronger in Fig. 4c for 7 = 1 than for 
7 ^ 1 in Fig. 4a. Similarly, a much stronger variation 
of Ic{T)/Ic{0) is seen in Fig. 4b for 7 = 1 than for the 
(not pictured) BCS-like results obtained with 7 <C 1. 
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FIG. 4. Plots of the normalized critical current 
Ic{T) / IciQ) for superconductors with Fermi surface FS2 and 
Jn ~ 1,25,50,100 meV. (a) d-wave and extended-s with 
7 <C 1. (b) s-wave with 7 = 1. (c) d-wave and extended-s 
with 7=1. 



We remark that Eq. (|8|) is a fully general expres- 
sion for the c-axis tunneling between two layered super- 
conductors, assuming that only one OP component ex- 
ists in each superconductor, and that all tunneling pro- 
cesses are completely coherent. That is, the L and R 
superconductors can have different quasiparticle disper- 
sions (both parallel and normal to the junction), different 
OP symmetries, and different Tc values. Thus, Eq. ( p^ ) 
also describes the tunneling between a three-dimensional, 
conventional superconductor placed on the top of a lay- 
ered superconductor with an unknown OP symmetry. It 
also applies to the case of coherent tunneling between 
two three-dimensional superconductors. When one of 
the superconductors is a three-dimensional material, one 
simply modifies the quasiparticle dispersion ^^(k, fc^) — 
^Ot]{^) — Jt] cos(fczs) far from the junction, by setting J,j| 
sufhciently large, and comparable to J^||. Except for 
the case of real-space pairing with a component normal 
to the junction, the fact that the pairing was assumed to 
take place only within the layers is not important, since 
the layer index is dropped, so that the pairing strength 
is constant throughout each superconductor. Thus, a su- 
perconductor with an isotropic, three-dimensional quasi- 
particle dispersion would have ^,,(k, fcj) = (k^-|-fc^)/(2m) 
far from the junction, where m is the quasiparticle ef- 
fective mass. This form is obtained from Eq. (2) by 
setting J^ii — (ma^)^^ and J,, = (ms^)^^, and requir- 
ing J^, J,j|| » ep- An anisotropic three-dimensional su- 
perconductor might have either or J^|| different from 
these values. 

We remark additionally that for the two-dimensional 
tight-binding Fermi surfaces, one restricts the wavevec- 
tor integration to the first Brillouin zone (BZ), which for 
a tetragonal material has Ifc^l, < n/a. In the three- 
dimensional limit, one can use the free-particle form for 
the quasiparticle dispersion afr from the junction, and 
remove the BZ limits on the components kx,ky of the 
wavevectors parallel to the junction. Note that one still 
does not integrate over fc^, the wavevector normal to 
the junction, which is not a good quantum number in 
the presence of the junction. Thus, even in the limit of 
Josephson tunneling between two conventional bulk su- 
perconductors, surface states are present, and affect the 
tunneling. Except for the tunneling Hamiltonian limit 
7 << 1, for which we had previously included the ef- 
fects of these surfaces states, our present calculation 
is thus the first one to correctly include the surface states 
in coherent Josephson tunneling between two supercon- 
ductors, whether conventional bulk, or layered. 

It is important to note that when the quasiparticle 
dispersions in the two superconductors are not identical, 
there is an intrinsic impedance mismatch between the two 
materials. Since in coherent tunneling, the wavevector 
parallel to the junction is preserved exactly, if for some 
wavevectors the Fermi surfaces in the two half spaces 
are not identical, it is not possible for the tunneling at 
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those places in the BZ to be elastic. That is, one can- 
not preserve both the momenta and the energy. The 
result of this impedance mismatch is that the amplitude 
of the coherent tunneling, and hence the critical current, 
is reduced from what it would be if this impedance mis- 
match did not occur. This impedance mismatch should 
actually occur in coherent tunneling between all inequiv- 
alent superconductors. For example, in tunneling be- 
tween Nb and Pb, the Fermi surfaces, which are both 
three-dimensional in nature, are not identical, and some 
amount of coherent tunneling would be suppressed by 
this effect. The amount of the suppression ought to de- 
pend strongly upon the particular crystal surfaces stud- 
ied at the junction location. 



VI. C-AXIS TWIST JUNCTIONS 

We now consider the case of purely coherent c-axis tun- 
neling between identical layered superconductors twisted 
an angle 4>o about the c-axis with respect to each other. 
This is a special case of coherent tunneling between two 
different layered superconductors. However, the OP sym- 
metry must be the same on each superconductor, with 
the only caveat that the OP has to arise from the local 
pairing interaction, Eq. ( 14 ) . The only difference is that 
the crystal orientation is rotated by ±0o/2 in the two 



half spaces, and this leads to similarly rotated quasipar- 
ticle wavevectors. Thus, in the two superconducting half 



spaces, we take k,. 



(fcj;^,fcy^), which are rotated by 



±</)o/2 for r] = L,U, respectively. For any OP symme- 
try, we thus have ^or;(k) = Co(k,,) and A^(k) = A(k^). 
Although for an isotropic s-wave pairing interaction, the 
twist orientation doesn't matter as far as the OP A is 
concerned, a non- vanishing (/>o still causes the quasiparti- 
cle dispersions to be different in the two superconductors. 
This leads to a strong impedance mismatch for all 00 val- 
ues not too close to 0°, 90°. 

In Figs. 5-8, we have presented our results for 
IciM/IciO) at T/Tc = 0.5 for different values of the 
material parameters. In Fig. 5, we presented our re- 
sults for the optimally doped tight-binding quasiparti- 
cle dispersion, Eq. (p^, in which /i = 0.6, v = 1.3, 
and J[[ = 500 meV, with Fermi surface FS2. We in- 
clude curves for each of the three OP symmetries, and 
for Jjj = 25, 50, 100 meV. For the extended-s-wave OP, 
we set e = 0, so that the wavevector dependence of the 
OP has the | cos(fc2,a) — cos(fcj^a)| form. Results for the 
extended-s-wave OP with e > are intermediate to the 
extended-s and s-wave results shown. 



100 "S-wave fn\ 

50 A d-wave 
25 X extended-s 




10 20 30 40 




1 20 30 40 

FIG. 5. Plots of 7c(</>o)//c(0) at T/T^ = 0.5 for c-axis twist 
junctions between layered superconductors with FS2 and OPs 
of S-, d-, and extended-s-wave symmetry, for Jrj — 25, 50, 100 
meV. (a) 7 <g 1, including Umklapp processes, (b) 7 <^ 1, 
without Umklapp processes, (c) 7 = 1 with Umklapp. 
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In each case, there are both direct and Umklapp tun- 
nehng processes. In the direct processes, a quasiparti- 
cle undergoing tunnehng across the twist junction has 
wavevectors and k^/ both within the first BZ. How- 
ever, for the case in which either the initial or the final 
wavevectors is outside of the respective first BZ, Umk- 
lapp processes can occur, due to the periodic nature of 
the quasiparticle dispersions assumed. To investigate 
whether these Umklapp processes are important, we treat 
the two limiting cases, either (1) that they can be com- 
pletely neglected, or (2) that they are equal in weight to 
the direct processes. 

In Fig. 5a, we set 7 ^ 1, the tunneling Hamiltonian 
limit. In this case, we included the Umklapp processes 
with the same weighting as for direct tunneling processes 
within the first BZ on each side of the twist junction. 
Preliminary versions of Fig. 5a with slightly different 
parameters were presented earlier. [|l0| For comparison, 
in Fig. 5b, we used the same parameters as in Fig. 5a, 
but the Umklapp processes were completely excluded. 
We note that for this quasiparticle dispersion, not very 
much of the intersection of the Fermi surfaces is excluded 
by neglecting the Umklapp processes for a twisted junc- 
tion. In addition, in Fig. 5c, we presented our results 
for 7=1, including the Umklapp processes. We note 
that comparing these results with those of Fig. 5b, the 
main differences occur for the larger values, especially 
Jn = 100 meV. Otherwise, for small J,,, there is very 
little difference between them. 

In addition, we note that the main differences between 
the (i^2„j^2-wave and extended- s- wave Ic{4'q) results ap- 
pear for 00 close to 45°. As 0o ^ 45°, Ic{4>o) I Ic{'^) 
for the d-wave OP, whereas for the extended-s-wave OP, 
Ici4'o) / lc{0) remains finite and flattens out, becoming 
only weakly dependent on ^q. Note that for Fermi surface 
FS2, even the isotropic s-wave OP leads to an anisotropic 
Ic{(j>o) / lc{0), for each of the Jjj values shown, reflect- 
ing the impedance mismatch across the twist junction. 
The greatest (/)o variation occurs for the smallest J,, val- 
ues, as the quasiparticle dispersions are the most two- 
dimensional, with the greatest impedance mismatch oc- 
curring for (po far from 0° , 90° . 

In Fig. 6, we again plotted /c(0o)//c(O) at T/Tc = 0.5 
fro c-axis twist junctions between optimally doped sani- 
pels with Fermi surface FS2 in the tunneling Hamiltonian 
limit 7 <^ 1, including the Umklapp processes. But, we 
now used the wavevector-dependent coherent interlayer 
tunneling similar to that suggested by Liechtenstein et 
al. 10] We note that this wavevector-dependent interlayer 
tunneling gives rise to a stronger 0o dependence of Ic 
than does the constant, wavevector-independent tunnel- 
ing model. In addition, for this model, the cxtendcd-s 
and d-wave OPs give rise to nearly identical Ic{4'o) / lc{0) 
curves, except for (jjQ near to 45°, of course, where 
Ic{4>o) for the d-wave case, but not for the other 
two OPs. 



0.8 



° 0.6 



S 0.4 



0.2 



•, \ X J(k,k')= 


FS2 

J(p(k)(p(k') 




s-wave 

extended-s 

d-wave 


I ^ 


200 
50 

\ pc 






"■---1::;:-'--'* 





10 



20 



30 



40 



<t>n (deg) 



FIG. 6. Plots of IciM/IciO) at T/Tc = 0.5 and 7 < 1 
with Umklapp for c-axis twist junctions between layered 
superconductors with FS2 and J(k, k') = Jip(k)ifi{\i), 
J,,(k,,,k^) — Jr,(p^(kr,), where v'(k) = \ cos{kxa) — cos{kya)\, 
kj7 = k and kz, — k', = 25,50,200 meV and s-, d-, and 
extended-s-wave OP symmetry. 

In Fig. 7, we plotted /c(0o)/^c(O) at T/Tc = 0.5 in the 
tunneling Hamiltonian limit (7^1) for c-axis twist junc- 
tions between heavily underdoped samples with Fermi 
surface FSl, including Umklapp processes, and kept 
fixed at 100 meV. Although curves are shown for each 
of the three OP symmetries studied, for this quasiparti- 
cle dispersion, it is difficult to distinguish them, as the 
Ic{4'o) / lc{0) values all either vanish (for the d-wave OP) 
or become very small for (j)o ~ 45°. 
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FIG. 7. Plots of /c((?!)o)//c(0) at T/Tc = 0.5 and 7 < 1 with 
Umklapp for c-axis twist junctions of layered superconductors 
with FSl with Jr, = 100 meV, and s-, d-, and extended-s-wave 
OP symmetry. 
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Finally, in Fig. 8, we presented plots of Ic{(j)o) / lc{0) at 
T/Tc = 0.5 in the tunneling Hamiltonian limit (7^1) 
for c-axis twist junctions between heavily overdoped sam- 
ples with Fermi surface FS3. Results for = 25, 50, 100 
meV are shown. In this case, there is very little 
impedance mismatch for the isotropic s-wave OP, since 
the rotation induced by the twist junction nearly maps 
the Fermi surface onto itself. One of these curves (the s- 
wave OP curve with J,, = 100 mcV) is actually consistent 
with the data of Li et ai, although the quasiparticle dis- 
persion assumed is very metallic and three-dimensional, 
completely inconsistent with that expected for BSCCO. 
The s-wave curve with = 50 meV is only marginally 
consistent with the data at best, and the s-wave curve 
with smaller values can be ruled out. Similarly, the 
extended-s and especially the d-wave curves are all in- 
consistent with the experimental data. 




FIG. 8. Plots of h{(l>o) / IciO) at T/Tc = 0.5 and 7 < 1 with 
Umklapp for c-axis twist junctions of layered superconductors 
with FS3 and with s-, d-, and extended- s-wave OP symmetry. 

VII. CONCLUSIONS 

We calculated the Josephson current for c-axis coher- 
ent tunneling between two layered superconductors ex- 
actly. We assumed the intralayer quasiparticle dispersion 
of the tight-binding form, Eq. Since the parame- 

ters in this model are arbitrary, our results also apply to 
the cases in which one or both of the superconductors is 
a conventional, bulk material. In these cases, one sim- 
ply lets one or both of the ^ Jj^W — * 00. For 7 7^ 1, 
the tunneling properties across the junction are different 
than the intrinsic tunneling between adjacent layers far 
from the junction, and the role of surface states on the 
sides of the junction becomes important. 

Such surface states were first investigated for super- 



conductors in the case of a vacuum interface. Then, 
the role of surface states upon the coherent Josephson 
tunneling between layered superconductors was investi- 
gated within the tunneling Hamiltonian limit. |^ Our 
present results comprise the first calculation of the strong 
coherent Josephson tunneling between two superconduc- 
tors which properly takes account of the surfaces states. 
The earlier results of Tanaka and Kashiwaya, 0| which 
were claimed to give the correct results for this prob- 
lem, are completely incorrect. In those calculations, the 
quasiparticle dispersions were not of the proper layered, 
tight-binding form, and the approximations used were 
not correct in any limit. 

Our results also apply to the case of tunneling between 
two superconductors in other orientations. For instance, 
in the case of an isotropic bulk superconductor, the sur- 
face states present along the sides of a junction normal to 
the c-axis are precisely the same as that along the sides 
of a junction in another orientation. Thus, the calcula- 
tions done by others, including Tanaka and Kashiwaya, 
for the Josephson tunneling within the a6-plane of 
layered superconductors, do not correctly take these sur- 
face states into account. In addition, since they always 
assume a circular Fermi surface cross-section (similar to 
FS3 in Fig. 2), all effects of impedance mismatching 
present in our investigation of the c-axis twist junctions 
are neglected. As shown in Section VI, such effects can be 
very strong, especially when the OP is highly anisotropic. 

Although we have not introduced the effects of inco- 
herent c-axis tunneling, we found that coherent tunnel- 
ing between two superconductors that are not precisely 
the same gives rise to strong impedance mismatch ef- 
fects. For coherent tunneling between identical layered 
superconductors twisted about the c-axis with respect to 
each other, the impedance mismatch effects are strong. 
Even when the OP is isotropic, one expects a strong twist 
angle (po dependence of the Josephson critical current. 
Although Umklapp processes are present in tunneling 
between superconductors twisted about the c-axis, they 
make only a very small correction to the actual critical 
current. 

We calculated the Josephson critical current across c- 
axis twist junctions, for OPs of the isotropic s-wave, the 
tight-binding d^2_y2--wave, cos{kx:a) — cos(fcj,a), and the 
tight-binding 'extended-s-wave', {[cos(/c2;a)— cos(/cya)]^-|- 
g2|i/2^ forms. We studied the cases of relative weak tun- 
neling (7 <C 1) and strong tunneling (7 = 1), for various 
c-axis dispersion bandwidths. Except for the simple case 
of an isotropic OP on a circular Fermi surface (akin to 
FS3), we conclude that coherent tunneling cannot possi- 
bly explain the data of Li et al. ^ Moreover, since the 
dx2-y2-'wsLve OP cannot possibly fit the data of Li et al. 
under any circumstances, we conclude strongly that the 
OP is not c?-wave near to Tc- Since there is widespread 
agreement that the Fermi surface of BSCCO is unlikely 
to have a circular cross-section, our calculations thus sup- 
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port the contention of Li et al. that the c-axis tunnel- 
ing intrinsic in BSCCO must be largely incoherent. ||] 
A brief discussion of the role of incoherent tunneling has 
been presented, supporting this argument. | |lO| A more 
thorough treatment of the role of incoherent tunneling 
will be discussed in detail elsewhere. llTlll 
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cos(6',,±) 



(26) 



The parameter 9^+ is defined to have a positive imag- 
inary part, and Orj- is defined to have negative imag- 
inary part. For convenience, we have defined f2,j = 
^uP- ^ I Pi Ci; = and b„ = A^/fi^. For com- 

parison with previous notation, M we then write 



exp(i6l^+) = yS^, 



(27) 



where is given by Eq. (||). 

We then find the half-space Green's functions, using 
the procedure of Eqs. (|^) and (^). For the L half space 
defined by m, n < 0, we have 



IX. APPENDIX: EXACT SOLUTIONS FOR 
GREEN'S FUNCTIONS IN THE LAYERS 

To obtain the elements C/,„„ of the bulk tight-binding 
Green's function matrix we can first Fourier series 
transform in the layer index, letting 



sdk. 



TT /s 



2tt 



Eq. (|l|) can then be solved algebraically to give 

iuj + ArjTi + ^r/Ts 



e''(k,fc„c^) 



£2 



where the quasiparticle dispersion 

£,,(k, fc^) = Cor,(k) - Jjjcosk^s. 



(19) 



(20) 



(21) 



For notational simplicity in the following, we shall omit 
the k and lu dependencies of the various functions. 

Performing the integral in Eq. (|l9|) and multiplying 
by {Jri/2)T3 on the right, we write 



L — qL _ qL 'li^ L 

Ivan =^mn ^ml r) ^iJGn^ 



(28) 



which is easily solved for g^^ by setting to = 0, yielding 

gLu = Qtn - + Qoi'^^^VqL (29) 

Similarly, for the U half space defined by TO,n > 1, we 
find 

9mn = Gmn ^ iGmO-^'^3)[^ + GlQ-^Ts]^^ ^in- (30) 

After some straightforward algebra, one finds for each 
half space 



where 



and 



^ hi + n't . f 
— ^mn ' mil ' 



(31) 



(32) 



where 



and 



with 



Pin = (e+„ + en)/2, 



rnn 



±j exp(±i0,,± |n — to|) 
2sin(0^±) 



(22) 

(23) 
(24) 

(25) 



^mn 2 VAmn ' AmnJ 

for rj = L,U. For n < to < 0, we have 

L± _ - cxp[TiOL±{n - 1)] sin[6>L±(TO - 1)] 

For 1 < TO < n, 

exp(±i0(7±n) sin(0(7±TO) 



sm{eu±) 



(33) 



(34) 



(35) 



To allow for different OP phases in the two sides of 
the junction, we introduce a phase factor multiplying the 
order parameter in the left hand side, letting 



Sl = exp[i(0L - M]- 



(36) 
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The function Su thus does not inchide such a phase fac- 
tor. 

So far, we have found expressions for the two distinct 
half space Green's functions, which are electronically un- 
coupled from each other, since no quasiparticle propa- 
gation from one half space to the other has yet been 
introduced. We thus couple them together via the local 
perturbation J with matrix elements given by Eq. (|^). 

The exact solution to this problem of coupled half 
spaces then yields the full Green's function matrix G, 
with matrix elements Gmn- In the L half space, n, m < 0, 
the Gmn satisfy 



G. 



Jl 



^nn^Tz= 9mn+19tQ[{9ii) ^ " 75oo] ^ qL- (37) 



Now the gap equation, Eq. (0) in Section 111 can be 
found from G„„. On the L side of the junction, we then 
need Gqo- From Eq. ( p7| ) wc find that the exact Green's 
function at the interface and in the lower superconductor 



Combining, and evaluating the trace, we find 
\5l\\5u\AuAl 

k LO 



D 



sin(0£/ -(t>L). (45) 



From this expression, it is elementary to obtain Eq. ( p^ ) 
of the text. 

In the tunneling Hamiltonian limit, 7 << 1, this re- 
duces to 

/ = -4e7T^^ |,5i||(5c,|3[exp(*0t,+ )]x 



k 



x5[exp(i0i+)] sin( 



(46) 



This result agrees with the tunneling Hamiltonian result 
for the coherent tunneling limit, as expected. [H 



where 



00^0 - 7l exp(z0t;+ + i6L+)?{gYi)-' 



D 



(38) 



I? = l + 72|exp(i0t/+-Kz0L+)|' 

-2'y{[eL^u + h^L cos{(f>L ~ (j)u)]AuAi 
+BuBl}, 



and 



Ar, = [exp(i6l^+) - exp(-i6'^_)]/2 



Brj = [exp(i6l^+) -I- exp(-i6'^_)]/2. 



(39) 



(40) 



(41) 



Then, the trace in Eq. (|^) can be evaluated by first 
multiplying Eq. ( ^ ) by 2t^/ Jl on the right side. Leading 
to 



Tr[(Ti+ir2)Goo] = ^[AL5L 



|exp(i6'i+)pAc;(5t 



(42) 



This then leads to Eq. (9) in the text. 

In order to calculate the Josephson tunneling current, 
we need to find Gqi and Gio explicitly. Solving for these 
functions yields the exact results 



(43) 
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J 



Gio-Ta = 7[(gii) -75oo] 5oo- 
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